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Abstract: The classical two-body system with Lorentz- invariant Coulomb action- 
at-a-distance V = —k/p is solved in 3+1 dimensions using the manifestly covari- 
ant Hamiltonian mechanics of Stiickelberg. Particular solutions for the reduced 
motion are obtained which correspond to bound attractive, unbound attrac- 
tive, and repulsive scattering motion. A lack of perihelion precession is found 
in the bound attractive orbit, and the semi-classical hydrogen spectrum subse- 
quently contains no fine structure corrections. It is argued that this prediction 
is indicative of the correct classical special relativistic two-body theory. 



1 Introduction 

One of the most significant applications of the Stiickelberg relativistic dynamics 
[|], H is certainly the classical two-body problem. The significance arises from the 
fact that the Stiickelberg mechanics allows the use of Lorentz-invariant action- 
at-a-distance in the form V ~ V{p), p = p(t)|^, to model long-range mutual 
interaction between particles. 

In this paper we study the classical covariant Kepler problem, i.e., the clas- 
sical special relativistic two-body problem with scalar Coulomb action-at-a- 
distance V — —k/p, where A: > is an invariant constant. Although certain 
solutions for this potential were obtained by Cook |^ and by Piron and Reuse 
[|| , the full range of physical motion produced by this potential has not yet been 
investigated. 

In this paper we outline^ two separate methods to study the reduced two- 
body motion for this potential using the pseudospherical coordinates (p, 9, cf)) 

^The parameter r is the invariant dynamical time; see ref. j^. 

more complete discussion of these solutions is given in ref. [Q. The 1+1-dimensional 
version of this system is also discussed in ref. |fj] . 



suggested by CookQ Both methods are generahzations of solutions of the non- 
relativistic inverse-square potential V = —k/r. In the first method, we integrate 
the Euler-Lagrange equations for the central force V = V{p) to obtain a set of 
Lorentz-invariant isolating integrals, which are solved for a set of classical orbit 
equations in the pseudospherical coordinates. We plot particular solutions for 
the invariant Coulomb potential in the 2+1 dimensional particle separation co- 
ordinates (Aa;(T), Ay(r), At(T)). In the second method, we solve the classical 
Hamilton-Jacobi equation to obtain the Lorentz-invariant action variables for 
the central force potential. In the case of the Coulomb potential, we impose the 
rule of semi-classical quantization to obtain a prediction for the electron energy 
levels of hydrogen. 

The classical orbit solutions of the potential V — ~k/p will be shown to 
produce all three types of physical motion expected for the two-body system 
with inverse-square interaction — i.e., bound attractive, unbound attractive, and 
repulsive scattering motion — with the distinction being given in a straightfor- 
ward manner by a pair of Lorentz-invariant eccentricity constants^. The most 
striking prediction in these solutions is the lack of perihelion precession in the 
bound orbits of the special relativistic system^. As a result of this lack of 
precession, the semi-classical hydrogen spectrum contains no fine structure cor- 
rections. It is argued that this result is favorable to the Stiickelberg theory, 
based on consideration of the classical limit without spin. 



2 The Classical and Semi- Classical Solutions 

Consider an isolated system of two particles with constant finite rest masses m^, 
i — 1,2, and mutual long-range interaction over flat spacetime. It is well known 
[H ^ ^, |l^ that in the covariant mechanics, the most straightforward solution 
of the classical problem with scalar interaction is obtained by a transformation 
from the individual particle event coordinates {xi(t),X2{t)), fi = 0,1,2,3, to 
the covariant center- of-mass and reduced coordinates, {X^{t),x^{t)). The sys- 
tem can be regarded as comprising two fictitious particles with masses M and 
m. The transformations are 

AI — mi + m2, m — (mi 7712) /M, (1) 
= (toix^' + max^) /M, x^" = x!^ - (2) 

In these coordinates, the nontrivial motion of the system occurs entirely in the 
reduced problem, for which the frame components are x^ — {At, Ax, Ay, Az). 

The most convenient procedure is to make a second coordinate transforma- 
tion of the reduced problem to the pseudospherical coordinates of Cook 
which consist of the variable spacetime separation p = p{t), as well as a space- 
time angle /3 = /3(t), and two ordinary angles 9 = 9{t) and (j) = 4'{t)- The 

•^In ref. 1^1; see Section 3 below for a discussion of these coordinates. 

*Here we set A: > and obtain both attractive and repulsive solutions. The distinction 
between the two cases is absorbed in the definition of the eccentricity constants; see ref. |p|. 
5cf. A. Sommerfeld, 1915 ^ 
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transformations are 



(Ax)^ + {^yf + (Az)2 



1/2 



tanh/? M/d(l\xY + {l\yY + {l\zY 



(3) 



tane* = Az/VAa?2 + A?^ 
tan</) = Ay/Aa;. 

In these coordinates, the Lagrangian for a central force potential is 



0^ (4) 



where the dot indicates differentiation with respect to the world time. By in- 
tegrating the four resultant Euler-Lagrange equations in the pseudospherical 
coordinates {p, (3,9,(j)), one obtains the Lorentz-scalar isolating integrals, 



lcl> 

k 

E 



mp cosh/30, 

p4 cosh^ f30^ + ll csc^ 

mV"/?^ - ll sech^/?, 
, A2 



1/2 



(5) 



1 .2 

-mp 



2rap-' 



V{p). 



It can be shown[^ that the conservation of a covariant angular momentum 
tensor implies that the polar angle 6 can be eliminated from the problem by 
an ordinary three- rotation of the coordinate axes, i.e., by choice of a frame in 
which = and 6 = 7r/2. Since this rotation does not involve relativity, it 
then possible to study the classical problem in 2+1 dimensions without loss of 
generality. 

By comparing differential expressions, the three remaining integrals of the 
motion (Z^, A^, E) may be solved to derive a pair of central force orbit equations 

for the dependence of p and (j) upon the spacetime angle /3. The azimuthal 
orbital component (p = <j){P) is given by 



m 



tan 




(7sinh(/3 - fip) 
-cosh2(/3-/3o) 



(6) 



where = — Z^/A^. For a conservative potential V = V{p), the radial orbital 
component p = p{/3) is given by the solution of the differential equation, 



cosh^/3/ d(3' 



, tanh f3 du 
cosh^ (3 'dp 



m dV 
A2 d^' 



(7) 



where u= 1/ p. 
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In the case of the two-body Coulomb potential, 



— k 

V — , k = Lorentz-invariant constant, (8) 

P 

the inhomogeneous term in the radial orbit equation is a constant, and 
thus particular solutions are known trivially once the homogeneous equation is 
solved. It is useful to examine the two particular solutions formed by adding 
the inhomogeneous term to the two linearly independent solutions of the ho- 
mogeneous equation, which may be labeled Type I and Type II solutions for 
P = p{P), 

Type I: z = ^(l-Tsinh(/3-/3o)), (9) 
Type II: 



1 


mk 


(-7 


p 




1 
p 


mk 


(-7 



(sinh2(/3-/?o)-/2)'/']) , (10) 



where / = -y/g^ — 1, and where e' are e" are two new invariant constants which 
are explicit functions of E and A Q . They may be regarded as Lorentz-invariant 
generalizations of the Kepler eccentricity [|ll| . Note that / is the same constant 
in both (|) and 

From the orbital solutions 9 = 9{(3) and p = p{P), the physical motion 
produced by the Coulomb interaction can be understood by using the inverse 
transformations of (^) to plot in the frame coordinates (Aa;(/3), Ay(/3), At(/3)). 
If it is desired, the r dependence of the solution may be included by numerical 
inversion of the expression r = r(/3) obtained from integration of eqs. (H) 

It is useful to examine the dependence of the Type I and Type II solutions 
upon the eccentricities e' and e", as well as upon /, the rotational constant]^ 
For a given e' or e", the effect of the variation of / is simply to modify the 
scale and proportion of the curve (see j^). The topology of the reduced orbit — 
and thus the classification of the physical motion — is found to depend on the 
magnitude of e' or e". 

Fig. (|l|) is a summary of the classification of orbits arising from variation of e' 
and e". For a proof of the properties discussed here, see ref. Q. From the three 
panels of the diagram, it may be seen that the Type I solutions correspond to 
an attractive Coulomb interaction, whereas Type II corresponds to a repulsive 
potential. For the Type I solutions, the effect of the transformation e' — e' is a 
reflection of the curve across the (Ax, Ay) plane, and thus one need consider only 
|e'|. For the range < e' < 1, the motion is bound, and thus a generalization 
of the elliptical Kepler solutions. For e' > 1, the motion is unbound, i.e., a 
generalization of hyperbolic Kepler solutions. The boundary case e' = 1 is 
bound at p ^ oo. The degenerate case e' — e" = corresponds to the limit 
TO ^ 0, M < oo For the Type II solutions, eccentricities in the range e" > 1 
produce a pair of conjugate curves; the reduced trajectory follows one of the 
two branches. Solutions in the range e" < 1 are bound solutions which are not 

^So-called because it depends explicitly on in 2+1 dimensions. 
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difFerentiable at </> = (0, tt) and thus may be considered unphysical. Likewise the 
general solution e', e" ^ appears to be unphysical on the same grounds based 
on numerical investigation |^. The Type I solutions were found by Cook Q 
who discussed only the bound case. The bound Type I orbits were also derived 
by Piron and Reuse Q| by a transformation out of the center-of-mass rest frame 
to one in which At(T) = 0. 

The most interesting feature of the solutions is the fact that the perihelion 
of the bound two-body system does not precess. This follows from the fact that 
the bound orbit of any central potential is a closed curve M . This is contrary to 
the prediction of a finite precession made by Sommerfeld |[7|] using the Heaviside- 
Lorentz equations in the one-body limit, i.e., infinite mass of the source. This 
latter precession is entirely special relativistic, and in the case of gravitation, the 
Sommerfeld rate is exactly one-sixth the rate predicted by general relativity for 
the same planetary masses. Sommerfeld proved that the nonvanishing perihe- 
lion precession results in the "correct" fine structure hydrogen spectrum under 
the rule of semi-classical quantization 0. This result has become questionable, 
however, and it is generally regarded today as incorrect This judgment is 
based primarily on the fact that in the full quantum theory, fine structure is 
known to arise intimately from the spin of the electron p3[ |, which has no clas- 
sical counterpart. One should expect, therefore, that in the classical limit, the 
special relativistic two-body system should not undergo perihelion precession. 

It is interesting to impose the rule of semi-classical quantization on the bound 
Type I solutions above. In the case of the hydrogen atom, one should expect the 
result to be the Bohr prediction without fine structure corrections. Since the 
energy spectrum is a function of the separation of energy states, it is possible 
to continue to ignore the energy of the center-of-mass and to proceed by the 
reduced problem alone. 

First it is necessary to solve the covariant Hamilton-Jacobi equation to obtain 
the invariant action variables. As in the Lagrangian version above, many of the 
results are valid generally for the conservative potential V — V{p). Following 
Goldstein ch. 10, it is useful to set aside for the moment the knowledge 
that the system can be reduced in dimensionality. Using the 3-|-l-dimensional 
Lagrangian in the canonically conjugate momenta are 

Pp = mp, pp ^ -mpP, 

Pe = mp cosh 130, p^ — mp cosh (3 s\n.04>- (11) 



The Lorentz-invariant Hamiltonian is 

1 



2m 




cosh^ (3 



A. 

SVC? t 



V{p). (12) 
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Using the rule pp — + dS/dp, etc., the energy equation K — E becomes the 
Hamilton- J acohi equation^ 



dp 



dp 



1 



dSe 
89 



cosh /? 

1 ros. 



sm 



2mV{p) = 2mE. (13) 



Hamilton's characteristic function S is assumed here to be totally separable, 
i.e., in the form S = Sp + Sp + Sg + S^. This allows one to solve eq. ( |T^ ) by a 
procedure which is a generalization of the nonrelativistic method for the central 
force problem Jill , i.e., using the invariant separation constants, (a^, ag,ap, E) 
M . This procedure yields the solutions 



dS 



"0, 



dS_ 

as 

as 

dp 




(14) 



cosh^/? ^' 



' 2mi; - 2mV - 



The action variables are given hy Jp = § dS/dpdp, etc., where the integration 
is taken over the full range of oscillation of the coordinate. The first two integra- 
tions are solved entirely in the manner of the nonrelativistic problem, yielding 
= 27ra0 and Jg = 2'k (ag — a^) The new step is the third integration (see 
[p|), which yields Jp — 2t: {ap — ag). The last integration is over libration of the 
radial coordinate p. The form of the integral is identical to the corresponding 
nonrelativistic one for r, and thus the identical method of contour integration 
is used ll^, yielding 



Jp = - (Jp + J9 + J4>) + 



2'Kimk 



V2mE 

It follows immediately that the energy is given by 



K 



[Jp + Jf3 + Jg + J^Y 



E. 



(15) 



(16) 



The angular frequencies, which are given by the rule Vp = dK/dJp, etc., are 
therefore identical. Thus it is proved what was asserted above, namely that the 
bound orbits of the central force potential are closed. 
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The rule of semi-classical quantization stipulates that Jp = Uph, etc., where 
the quantum numbers {np,np,ng,n^) take on integer values and where h is 
Planck's constant. From eq. (p^, the semi-classical energy levels in the case of 
the potential V ^ —k/p are 

^ = -^^17^^' «-l,2,3,..., (17) 
n = Up + np + Ufjj + ne- (18) 

Assuming that k represents the electrostatic interaction strength between the 
proton and electron, the prediction of the Stiickelberg special relativistic two- 
body theory gives the Bohr levels without fine structure corrections, which is 
the correct result in the case of no spin. 



3 Comments on the Coordinate System 

The pseudospherical coordinates eqs. (^) used in this paper are those used by 
Cook 1^ to study several two-body systems, including the inverse-square poten- 
tial V = —k/ p. These coordinates may be contrasted with the alternate set of 
pseudospherical coordinates used by Arshansky and Horwitz iQ to solve the 
full quantum Stiickelberg model of the hydrogen atom with no spin. Although 
Cook claimed to have found in the full quantum case an incorrect n + 1/2 de- 
pendence of the hydrogen spectrum on the principal quantum number n, this 
result is obtained by using an electromagnetic retarded potential approximated 
for low particle acceleration (cf. eq. (81) on p. 133 of ref. (1^). The inverse- 
square potential V = —kjp studied here is discussed by Cook only in the purely 
classical context, as a possible special relativistic gravitational solution. He 
used the bound two-body solutions he obtained for this potential to study high- 
order deviations from Newtonian motion in the Solar System orbits. The lack 
of perihelion precession discussed here is implicit in his work. 

Regarding the lack of fine structure corrections for the system without spin, 
the semi-classical prediction obtained in this paper is in agreement with the 
full quantum Stiickelberg result for this potential obtained by Arshansky and 
Horwitz using the alternate pseudospherical coordinates. In the purely classical 
case, the difference between the two coordinate systems is trivial, since the two 
sets are identical in 2-1-1 dimensions. In the quantum case, it may be possible to 
distinguish between the two coordinate systems based on the support required 
for the reduced wave function (see ref. |]l^). 
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Figure 1: Note for the electronic preprint version: this figure is viewable on- 
line as a GIF image at the address http://order.ph.utexas.edu/mtrump/figurell 
The classical orbits of the two-body problem for Coulomb action-at-a-distance 
V = —k/p, k > 0. The orbit shows the separation between the two particle 
events in the center-of-mass rest frame, using the Type I and Type II solutions 
in eqs. (^) and ( p^ as well as eq. (H). The Type I orbits correspond to an 
attractive interaction. For < |e'| < 1, the system is bound, as in fig. (a). The 
bound orbit is a closed curve which is differentiable; the nondiffercntiability in 
the diagram is due to numerical approximation. For |e'| > 1, the system is 
unbound, as in fig. (b). The Type II solutions, for the case e" > 1, correspond 
to repulsive scattering. The Type II orbits for e" < 1 are unphysical. It appears 
moreover that the general solution e', e" / is unphysical, based on numerical 
investigation. The orbit also depends on the constant /, but variation of this 
constant, as well as variation of e', e" within a specified range above, results in 
a change of the overall shape of the orbit but not in the type of motion. [The 
values of the constants used in the examples are: (a) e' = 0.5, / = 1.0; (b) 
e' = 1.5, / = 2.0; (c) e" = 1.5, / = 2.0] 
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http://arXiv.org/ps/physics/9805032vl 



